Plasma heating in the very early phase of solar flares 
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ABSTRACT 

In this paper we analyze soft and hard X-ray emission of the 2002 September 20 
Ml. 8 GOES class solar flare observed by RHESSI and GOES satellites. In this flare 
event, soft X-ray emission precedes the onset of the main buUc hard X-ray emission 
by ~5 min. This suggests that an additional heating mechanism may be at work at 
the early beginning of the flare. However RHESSI spectra indicate presence of the 
non-thermal electrons also before impulsive phase. So, we assumed that a dominant 
energy transport mechanism during rise phase of solar flares is electron beam-driven 
evaporation. We used non-thermal electron beams derived from RHESSI spectra as 
the heating source in a hydrodynamic model of the analyzed flare. We showed that 
energy delivered by non-thermal electron beams is sufficient to heat the flare loop to 
temperatures in which it emits soft X-ray closely following the GOES 1-8 A light- 
curve. We also analyze the number of non-thermal electrons, the low energy cut-off", 
electron spectral indices and the changes of these parameters with time. 

Subject headings: Sun: flares — Sun: X-rays, gamma rays 



-3- 



1. Introduction 



It is commonly accepted that during the impulsive phase of the solar flare, non-thermal 
electron beams accelerated anywhere in the solar corona move along magnetic field lines to the 
chromosphere where they deposit their energy. Here, most non-thermal electrons lose their energy 
in Coulomb collisions while a small fraction (~ 10"^) of the electrons' energy is converted into 
hard X-rays (HXR) by the bremsstrahlung process. The heated chromospheric plasma evaporates 
and radiates over a wide spectral range from hard X-rays or gamma rays to radio emission. 
These processes, called electron beam-driven evaporation model is believed to be a dominant 
energy transport mechanism during solar flares. One of the most spectacular and well observed 
manifestation of the described processes is soft X-ray (SXR) emission of magnetic loops observed 

lard and s oft X- ray fluxes 



Neupert 



(Il968h . Since 



e.g. by Yohkoh and Hinode satellites. Such a scenario implies that the 
emitted by solar flares are generally related, as was first described by 
HXR and microwave emissions are produced by non-thermal electrons and SXR are the thermal 
emission of a hot plasma, the Neupert eff"ect confirms that non-thermal electrons deliver energy 
spent on plasma heating. Thus the HXR emission is directly related to the flux of the accelerated 
electrons whereas the SXR emission is related to the energy deposited by the non-thermal electron 
flux. 



This interpretation, however, leads to a number of problems and quest ions concernin g 



energy content and tirn e rel ations between SXR and HXR fluxes (see e.g 



Dennis 



Dennis & Zarro 



(Il993h . and 



McTierman et al. 



mm- 



mm . 



There are several papers that inves tigate temporal d e pend enci es between the 



SXR and HXR emissions. For example. 



Machado et al. 



( 1986 ) and 



Schmal et al. 



jeginn ings of 



(Il989h reported 

frequent strong SXR emission before the impulsive phase. Following the se authors, on averag e, 



the SXR emission precedes the onset of HXR emission by about 2 min. 



Veronig et al. 



(l2002h 



analyzed 503 solar flares observed simultaneously in HXR, SXR and Ha. In more than 90% of 
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the analyzed flares, an increase of SXR emission began prior to the impulsive phase. On average 
the SXR emission starts 3 min before the hard X-ray emission. 

An enhanced thermal emission preceding the onset of the hard X-rays may be indicative of 
a the rmal prehea t ing ph ase prior to the impulsive electron acceler ation. Cu r rent- s heet models 



(e.g., 



Hevvaerts 



Li et al. 



(11987b discussed the 



(|1977n ) of solar flares predict a preheating phase, 
preheating phase of solar flares triggered by new emerging magnetic flux. They proposed four 
different types of reconnections to explain the preheating as well as impulsive phases of flares. 

Ano ther a pproach to this issue is multi-thread hydrodynamic modeling of solar flares. 



Warren 



(120061) suggested that modeling a flare as a bunch of independently heated threads may 
simulate precedence of the SXR emission. The author successfully reproduced the temporal 
evolution of high temperature flare plasma in the Masuda flare of 1992 January 13. 

An alternative mechanism to el ectron beam-driven eva poration, namely conducted driven 



evaporation, was developed lately by 



Battaglia et al. 



(|2009|) . They studied in detail the pre- flare 



phase of four solar flares using imaging and spectroscopy from the RHESSI satellite. These 
authors explain the time evolution of the observed emission for all analyzed events as an effect of 
saturated heat flux. 

The tendency of the soft X-ray flux to appear before ha rd X-rays emi ssion can be attributed 



also to the sensitivity threshold of the hard X-ray detectors (IDennis 



19881). At the beginning of 



the flare, the energy flux may be below the detection threshold of HXR emission. 

In this pape r we show, using unprecedented high sensitivity of the RHESSI detectors 



( Lin et al. 



20021) and a numerical model of flare, that early SXR emission observed prior to 
impulsive phase could be fully explained without any ad-hoc assumptions (at least for analyzed 
event). All necessary energy to explain the soft emission could be derived from observed HXR 
spectra. In Section 2 we describe the analyzed event. Section 3 presents the details of the HXR 
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spectra fitting, numerical modeling and the results. The discussion and conclusions are presented 
in Section 5. 

2. Observations 

For our work we selected RHESSI observations obtained without the activation of the 
attenuators in order to prevent discontinuities in fitting parameters. A solar flare with a simple 
single-loop structure was then chosen for convenience in numerical modeling. The investigated 
flare occurred in the southwest hemisphere in active region NOAA 10126 (S23E69) on 2002 
September 20. It was classified as an Ml. 8 GOES class flare. GOES X-ray light-curves of 
the flare are shown in Figure 1 (top panel). GOES (1-8 A) flux has a background level of 
1.04 X 10"^ W/m^ (C1.04). The SXR emission started to increase slowly at 09:18:15 UT and 
showed two local maxima at 09:21:00 UT and 09:22:30 UT. It reached its maximum at 09:28:30 
UT. Harder (0.5-4 A) GOES emission started to increase at the same time as the softer one 
and peaked one minute earlier at 09:27:30 UT. It also showed two local, even more pronounced 
maxima of the emission. 

RHESSI X-ray light curves of the flare taken in four energy bands are shown in Figure 1 
(bottom panel). The impulsive phase in X-rays > 25 keV started at 09:25:24 UT and it had two 
maxima around 09:26 UT and 09:27 UT, respectively. In 25-50 keV energy range, a small spike 
of emission occurred between 09:24:16 UT and 09:24:32 UT. It appeared also as a small hump 
on the both GOES light-curves. The X-ray emission below 25 keV started to rise simultaneously 
with GOES emission. Three local peaks of emission are presented in these energy bands at the 
moments of SXR local maxima. 

Images obtained using RHESSI data with the CLEAN imaging algorithm revealed SXR 
emission at 6-12 keV and intermediate energy emission 12-25 keV coincident with the flare 
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20-Sep-2002 




09:16 09:20 09:24 09:28 09:32 09:36 
Time (UT) 



Fig. 1. — Top panel: GOES light curves in two energy bands : 1-8 A (upper line) and 0.5-4 A. 
The horizontal dotted line represents the preflare background level observed in soft channel. The 
arrow indicates the starting moment of the numerical model. Bottom panel: RHESSI light curves 
taken in four energy bands (from top to bottom): 4-12, 12-25, 25-50 and 100-300 keV. 

location. These observations indicate that SXR emission recorded by GOES in the early 
phase of flare came f rom analyzed event. Images obtained with the PIXON imaging algorithm 



(Hurford et al. 



3) 



20021) showed a single flare loop (see Figure 2). This image allows us to define 
geometrical proportions of the loop. The cross-sections of the loop S = 1.13 x lO^^cnP- was 
estimated as an area within a level equal to 30% of the maximum flux in the 25-35 keV energy 
range. Half-length of the loop Lq = 9.5 x lO^cm^ was estimated from the distances between the 
centers of gravity of the footpoints, assuming a semi-circular loop shape. 

Unfortunately SOHO/EIT telescope observed this AR before an after the flare. However 
images obtained with 195 A filter after event at 09:47:59 and 9:59:59 UT confirmed the single 
loop structure of the flare. 
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RHESSI 20-Sep-2002 09:26:42.000 UT 



1 1 1 II 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 II 1 1 1 ii|iii 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1| 1 1 1 1 1 1 1 1 1 1 1 



LDetectors: 2F, 3F, 4F, 5F, 6F, 8F, 9F 
; Energy range: 25 - 35 keV 
; Accumulation time: 8 sec 
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Fig. 2.— A 25-35 keV PIXON image with time accumulation from 09:26:42 UT to 09:26:50 UT, 
at maximum impulsive phase. An isocontour corresponds to 30% of the maximum flux. 



The RHESSI data were analyzed using RHESSI OSPEX package of the SolarSoftWare 
(SSW) system. Data were summed over the front segments of the seven detectors with detectors 
number 2 and 7 excluded. The spectra were measured with 4 sec temporal resolution in 158 
energy bands from 4 to 300 keV. We applied the energy widths dE = 0.3 keV within the range 
4-15 keV, dE= 1.0 keV in the range 15-100 keV, and dE = 5.0 keV above 100 keV. The analyzed 
spectra were corrected for pulse pile-up, decimation and albedo effects. We used full 2D detector 
response matrix to convert input photon fluxes to count rates. Before fitting spectra we removed 
the averaged non-flare background spectra. For energies below 50 keV, the background spectra 
were accumulated and averaged from pre-flare period between 09:00 and 09:06 UT. For energies 
above 50 keV we used a linear interpolations between the time intervals before and after the 
impulsive phase. 

In our model we used spectra taken after 09:18:15 UT when the SXR emission started to 
increase. We used 4 seconds-long time bins, but after background subtraction we increased the 



3. Calculations - modeling of the heating of the loop 
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accumulation time in the period from 09:18 UT to 09:23 UT to keep the positive counts rates in 
most of the energy bins in the 4-20 keV range. 

The spectra were fitted with single temperature thermal plus thick-target models (vth -l- 
thick). The thermal model was defined by single temperature and emission measure of the 
optically thin thermal plasma, and is based on the X-ray continuum an d line emission calculated 



by the CHIANTI atomic code (IDere et al 



are based on the coronal abundances of 



(1997) 



Landi et al. 



Feldman & Laming 



(120061) '). The element abundances 
(|200o|) . The thick target model was 



defined by the total integrated electron flux Nnth, the power-law index of the electron energy 
distribution 6, and the low energy cut-off of the electron distribution Figure 3 shows the fit 
of the spectrum, accumulated between 09:20:04 UT and 09:20:12 UT. In addition to the thermal 
component, there is an important power law shape for emission between 15 keV and 20 keV 
which can be recognized and fitted as thick-target emission of the non-thermal electrons. 

We assumed the electron beam-driven evaporation model of the solar flare. Therefore, we 
used in hydrodynamic model of the analyzed flare the non-thermal electrons beams derived from 
RHESSI spectra as the heating source via the Coulomb collisions. Heating was modeled using the 



approximation given by 



Fisher 



(119891) . In this work we used ID Naval Researc h Laboratory So lar 



Mariska et al. 



'lux Tube Model code by Mariska and his co-workers ( Mariska et al. 
1989D). This code was slightly modified by us. We included: new radiative loss and heating 



functions, the VAL-C model of the initial structure of the lower part of the loop, and double 
precision calculations. All the changes were teste d by comparison of the results with the results 



of the original model. For details see the paper by 



Falewicz et al. 



(|2009|) . An important problem 



which we meet during the modeling of the flares using original NRL code was an insufficient 
amount of the mater located in feet of the loops. To solve this problem we applied the VAL C 
model of the solar plasma, extended down using Solar Standard Model data. It was done solely 
in order to obtain big enough storage of the matter. All other aspects of the NRL model of the 
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Fig. 3. — RHESSI spectral fit results for data accumulated with a 8 sec time interval between 
09:20:04-09:20:12 UT. The spectrum was fitted with one temperature thermal model (blue color) 
and thick target model (green color) with energy cut-off" = 15.8 keV (please see the text for 
more details). Total fitted spectrum is shown in red color. The errors bars of the observed spectra 
are also shown. 

chromosphere are unchanged (radiation is suppressed, optically thick emission is not accounted 
for, and no account is taken of neutrals, the net conductivity flux is negligible). 

We modeled the evolution of the analyzed flare as follows: an initial, quasi- stationary 
pre-flare models of the flaring loop was built using geometrical and thermodynamic parameters 
estimated from RHESSI and GOES data. These initial parameters of the flaring loop were as 
follows: semi-length 9500 km, radius 1900 km, pressure in feet 22 dyn/cm^. Small volumetric 
heating was used to keep this model in the quasi stationary state on the pre-flare level of activity 
before the start of the non-thermal heating. Then we started to model the heating of the loop 
by non-thermal electrons adding a dose of energy and calculating the resulting GOES flux. We 
used thick target parameters A^„i/,, 6 and E^ obtained from fits of consecutive i^/ZfiSS/ spectra 
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for each time-step as input into the Fisher's heating function. However, an acceptable fit can be 
obtained for all Ec values in the range 5-30 keV for all following spectra. Because the low 
energy cut-off determines an amount of energy delivered to the loop, this non-uniqueness could 
be limited using an independent energetic condition, like observed 1-8 A GOES flux. Indeed, 
for each time step we adjusted Ec values in order to achieve conformity of the observed and 
modeled fluxes in GOES 1-8 A band. Such use of GOES 1-8 A flux put important limitation on 
the allowed low energy cut-off values and reduces importantly the non-uniqueness problem. 

Fig. 3 shows an example of the RHESSI fitted spectrum where the low energy cut-off 
Ec was adjusted to equalize synthesized and observed GOES fluxes in 1-8 A channel. The 
high energy part of this spectrum was finally fitted with the following thick target parameters: 
6 = 1 31 , Nnth = 7.2 X 10^^ electrons /sec, and Ec = 15.8 + 0.1 keV. This non-thermal electron 
beam contains total energy flux Fnth = 2.16±0.14xl0^^ erg/ sec and heated our loop during 8 
seconds between 09:20:04-09:20:12 UT. The heating increases the synthesized GOES class in 
this period from B1.93 to B2.11, which exactly corresponds to the changes of the observed (with 
background removed) GOES flux. Given errors of Ec and Fmh were derived by fitting the GOES 
1-8 A flux with an accuracy of 0.01 W/m2. Thermal component of the spectrum was fitted with 
emission measure EM = 4.25 X lO'^^cm"^ and temperature Te = 14AMK. Obtained temperature 
is consistent with temperature obtained from GOES data (ll.OM^. We fitted also the spectra 
with the purely thermal model, but obtained temperature of 20. 3M^ seems to be quite high with 
respect to the GOES one and makes us more confident in the interpretation with the presence of 
the non-thermal electrons well before the impulsive phase. 

The final result of our modeling is presented in Figure 4. The synthesized GOES 1-8 
A light-curve closely follows the observed one. The correspondence between observed and 
calculated fluxes in the case of the 0.5-4 A band is not so ideal. The inconsistency could be 
attributed to crude estimation of the initial loop conditions, errors in RHESSI spectra restoration. 
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Fig. 4. — Comparison of the observed (black thin lines) and calculated (thick blue line) GOES 
fluxes in 1-8 A (upper curves) and 0.5-4. A (lower curves) energy bands during initial phase of the 
analyzed flare. 

and the simplicity of our model. However, because the variations of the calculated 0.5-4 
A light-curve did not differ too much from the observed variations, it means that our model 
simulate the main physical processes in right way. As a result we fully restored the observed slow 
increase of the SXR flux recorded before impulsive heating using only non-thermal electrons 
characteristics derived from the observed HXR spectra. Thus we confirmed that in the analyzed 
flare variations of the SXR and HXR fluxes are consistent with the Neupert effect. 

Figure 5 presents time variations of the electron beam {thick target) parameters and RHESSI 
fluxes before the onset of the main bulk of hard X-ray emission at 09:25:24 UT. Non-thermal 
electron distributions, characterized by these parameters, delivered to the loop an amount of 
energy sufficient to produce the slow increase of the GOES emission observed before the 
impulsive phase. Low energy cut-off Ec changed during this period between 12 keV and 24 keV. 
Electron spectral index 6 varied between 4.5 and ^ 10, while total electron energy flux ranges 
from 3.3 x 10^^ ergs/sec to 1.4 x 10^^ ergs/sec. Obtained variations of the E^ reflect probably 
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Fig. 5. — Time evolution of the thick target parameters and RHESSI fluxes before the start of 
impulsive phase. From top to bottom: Nmh, 5, and observed fluxes at 12-25 and 25-50 keV. 

temporal variations of the processes in primary energy source region. On the other side obtained 
vari ations of the index 6 sh ow a clear general pattern of the soft-hard-soft spectral evolution (see 
e.g. iGrigis & Benz I (l2004h ). Additionally, all peaks in the 12-50 keV energy range are related to 
local increases in energy flux of the non-thermal electrons and so to local increases in heating of 
the loop. The two most pronounced peaks (and heating increases) at 09:20:16 UT and 09:22:16 
UT were manifested as local flux enhancements on both GOES light-curves. The maximum 
evaporation speed, achieved in a early phase of the simulation, was equal to 130 kinj s. The 
maximum temperature, density, and pressure obtained at the loop apex during the modeled period 
were equal to 13.5 MK, 1.3 x 10'' cnr^, and 462 dyn/cni^, respectively. 
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4. Conclusions 

The start and quick increase of HXR > 25 keV flux (defined as an impulsive phase of flare) is 
interpreted as an indication of the injection of the non-thermal electrons into the flaring loop and 
a beginning of the plasma heating by these electrons. However, often the SXR emission starts a 
few minutes earlier than the HXR which raises a question about the heating source in very early 
stage of the flare. Numerous solutions were proposed of which heating by thermal conduction 
and Multi-Stranded Loop model are the most well known. We showed that it is possible to fit 
soft (GOES) and hard (RHESSI) X-ray emissions of a solar flares well before beginning of the 
impulsive phase without any additional heating besides the heating by non-thermal electrons. This 
was made possible because of the unprecedented high sensitivity of the RHESSI detectors which 
are able to measure very low hard X-ray flux early in the flare. Part of the emission, mainly in 
the energy range < 25keV, is non-thermal in nature and indicates the presence of non-thermal 
electrons. In this case of a Ml. 8 GOES class flare, the non-thermal electron energy fluxes of the 
order of 10^^ ergs/ sec, derived under the thick target interpretation, fully explains the required 
heating of the plasma and resulting increase in SXR emission. 

The main limitations of our work are: crude estimation of the initial physical and geometrical 
parameters of the loop, errors in restoration of the RHESSI spectra and in GOES calibration, 
relative simplicity of the ID hydro-dynamical model and single loop approximation of the event. 
While the synthesized GOES 1-8 A light-curve closely follows the observed one, all limitations 
mentioned above caused small differences between observed and calculated fluxes in the 0.5-4 
A band. 

Various au t hors c onsidered classical heat conduction from the loop-top energy reservoir. 



Battaglia et al. 



(I2009|) fitted (RHESSI) spectra of four flares with a purely thermal emission and 
claim that at least for some flares electron beam heating does not work. In a case of flare analyzed 
by us purely thermal spectra seems to give too high temperatures (> 20 MK) comparing to those 
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obtained from GOES data. Additionally, even if the hard X-ray spectra observed with RHESSI 
provide no evidence for non-thermal p articles this does not necess arily means that the electrons 



don't exists. As it was pointed e.g. by 



Brosius & Holman 



(|2009|) the non-thermal hard X-ray 



emission associated with electron beam could be below RHESSI's level of detection. These 
authors estimated that in the case of analyzed event the non-thermal X-ray emission produced by 
an electron beam of sufficient energy flux to heat chromospheric plasma to the temperature and 
emission measure observed by RHESSI is below the observed background level. It is possible 
that for some flares heat conduction and for others beam driven evaporation works on the early 
phase of flare evolution. In our case the heat ing of the mat t er cau sed by non-thermal electrons is 



2-4 orders higher then local conductive flux. 



Stoiser et al. 



(|2008l) presented analytic predictions 



of the X-ray EM and expected in single and filamented flare loops for both mechanisms 
of evaporation and have tested these against real data consisting of 18 RHESSI microflares. 
Their results suggest beam heating in filamented loops to be in agreement with data. The peak 
temperatures were consistent with both single loop and multi-thread heating. The observed 
emission measures were mostly compatible with beam driving for a number of threads, but for 2 
events EM were also compatible with single loop model. 

Our results extend the standard model of SXR and HXR relationship to the early phases of 
solar flares and thus expands the number of flares consistent with the Neupert effect. These results 
also indicates that the process of electrons acceleration appear s during the early stage of the flare. 



Asai et al. 



(|2006h who reported 



well before the impulsive phase. This was confirmed e.g. by 
detection of the coronal non-thermal emission during the pre-impulsive phase of the X4.8 flare on 
2002 July 23. 

Our method of adjustment of the low energy cut-off" E^ in order to equalize synthesized and 
observed GOES fluxes in 1-8 A channel can be considered as a new method of Er determination. 
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